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Abstract 

We study the intermediate- and long-range forces between moving and spinning Skyrmions, employing 
two different methods. One uses a relativised product ansatz for the Skyrme fields, the other models 
Skyrmions as triplets of scalar dipoles. The methods lead to the same finite-dimensional Lagrangian 
dynamical system which may be interpreted as a point-particle approximation to Skyrmion dynamics. 
We discuss in detail the dynamics in the so-called attractive channel and the interaction between well- 
separated and rapidly spinning Skyrmions, and point out the resemblance between the latter and the 
one-pion exchange potential in nuclear physics. 



* e-mail: b.j.schroers@durham.ac.uk 



1 



1 Introduction 



The Skyrme model is a non-linear, non-integrable field theory, defined on 3+1 dimensional 
Minkowski space, which has soliton solutions 0. Here we will use the term "soliton" to mean 
what is also sometimes called a topological soliton: a localised, finite-energy and smooth solution 
of the field equations which is stable for topological reasons. Examples include lumps in the CP™ 
model, vortices in the abelian Higgs model and BPS monopoles in Yang-Mills-Higgs theory. In all 
cases the solitons have an associated integer winding number which is a topological invariant and 
which is also called the soliton number. Skyrme solitons of soliton number one are referred to as 
Skyrmions. Suitably quantised they are models for physical baryons. 

Given the long-standing difficulties in finding a satisfactory theoretical model for the interaction 
of two baryons (the nuclear two-body problem) it is not surprising that much effort has been 
devoted to the study of the classical and quantised interaction of two Skyrmions. In the absence 
of integrability one cannot hope to find exact time-dependent solutions which represent two- 
Skyrmion scattering or bound states, and one therefore has to resort either to numerical methods 
or to some kind of approximation. 

One scheme which has been used to great effect in the study of soliton dynamics at low energy 
is the adiabatic or moduli space approximation, explained systematically in ref. ||. In this 
scheme one views the field theory as an infinite-dimensional Lagrangian dynamical system (C, L), 
consisting of a configuration space C, the space of all fields at a given time, and a Lagrangian 
of the form L = T — V , where T is the kinetic energy functional and V is the potential energy 
functional. One then truncates the theory to a finite-dimensional Lagrangian system by choosing a 
finite-dimensional submanifold M. G C and restricting LtoAi. The space Ai is called the moduli 
space and should contain all configurations which, on general physical grounds, are deemed to be 
relevant at low energy. In theories with topological solitons the procedure is somewhat simplified 
because C decomposes naturally into disjoint sectors C n which are each labelled by the winding 
number n of the fields in it. Correspondingly, one may consider each winding number separately 
and choose, for given n, a truncation M. n of C n . The elements of Ai n are called n-solitons. 
The basic assumption of the moduli space approximation is that there exist finite-dimensional 
manifolds M. n such that the solutions of the truncated Lagrangian systems (A4 n , L\ Mn ) are, in 
some sense, good approximations to solutions of the full field theory. 

How does one find candidates for A4 n l M.\ is usually defined as the set of all fields with 
winding number one which minimise V. Typically, this set includes internal as well as position 
parameters. For n > 2, M. n should include configurations consisting of n well-separated single 
solitons with arbitrary internal parameters and positions, and one therefore requires that dim 
Ai n = nx dim M. x . In certain theories, which are said to be of Bogomol'nyi type, this requirement 
is automatically satisfied if one defines Ai n , like as the set of fields in C n which minimise V. It 
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follows that V is constant on A4 n , which means physically that there are no static forces between 
solitons of Bogomol'nyi type. As a result the truncated dynamics is entirely determined by the 
induced kinetic energy. This is generally of the form T = \gijC[iqj, where {g^} are coordinates on 
A4 n and the matrix g^, which depends on q i: defines a Riemannian metric on the moduli space. 
Thus in the moduli space approximation to the slow motion of solitons of Bogomol'nyi type the 
interactive dynamics of n solitons is modelled by geodesic motion on the moduli space. This is an 
adiabatic approximation because it approximates the time evolution of n-solitons by a sequence 
of static solutions of the field equations. So far the most interesting and successful application of 
moduli space approximation has been to two-soliton dynamics in theories of Bogomol'nyi type. 

Skyrme's theory is not of Bogomol'nyi type, but the static forces between Skyrme solitons are 
quite weak. It has therefore been suggested |2| that Skyrme solitons, like Bogomol'nyi type solitons, 
might be amenable to the moduli space approximation. This requires that one finds suitable 
moduli spaces for Skyrme's theory. A4i, defined as above, is six- dimensional and contains three 
translational and three rotational degrees of freedom. Thus AA.2 should be 12-dimensional. The 
set of minima of the potential energy V amongst two-Skyrmions is an eight-dimensional manifold 
of coincident Skyrmions || and does not include well-separated Skyrmions. On the other hand, 
there is a 12-dimensional family of two-Skyrmions, namely the fields of the product form, first 
considered by Skyrme in ref. |4j], which includes well-separated Skyrmions but does not correctly 
describe the low-energy coincident configurations. In ref. || N.S. Manton suggested a method 
of constructing a 12-dimensional manifold A4 12 which would include both the coincident and the 
well-separated Skyrme fields just described. More recently Atiyah and Manton have shown that 
one can obtain good approximations to Skyrme fields of low energy by calculating the holonomy 
of SU(2) instantons. They obtain a 12-dimensional manifold A4 12 of two-Skyrmions in this way 
and describe its topological and differentiable structure in ref. ||. They argue further that A4 12 
should be a good approximation to M 12 . The next step in the moduli space approximation, which 
has not yet been carried out, is to calculate the induced Lagrangian on A4 12 or on A4 12 . There 
are partial results for the potential, but little is known about the kinetic energy. 

In the present paper we study the kinetic energy of two well-separated Skyrmions which, 
statically, are well described by the product ansatz. Physically, this means studying the velocity- 
dependent forces between moving and spinning Skyrmions. According to the moduli space ap- 
proximation one should, in order to calculate these forces, allow the 12 parameters in the product 
ansatz to vary with time and express the kinetic energy in terms of the 12 parameters and their 
derivatives with respect to time. A few of the terms generated in this way were calculated in ref. 
H and a more complete numerical study was undertaken recently f7|. In this paper we will argue 
that this approach does not lead to the correct effective Lagrangian. Rather, one must allow the 
individual Skyrmions to Lorentz-contract when they move and modify the product ansatz to take 
into account this and other relativistic effects. 
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To back up our results we re-derive them using a different method. Our second calculation 
is based on the observation that, from afar, a Skyrmion looks like a triplet of point-like scalar 
dipoles. It has been known for a long time that the potential between two static Skyrmions 
can be understood in this way. Using retarded potentials we calculate the Lagrangian for the 
interaction of two spinning and moving triplets of scalar dipoles including all terms which are 
at most quadratic in linear and angular velocities. Similar retarded potential calculations, which 
may be interpreted as point-particle approximations to soliton dynamics, have been carried out 
for BPS monopoles |§ and maximally charged black holes || and have confirmed the moduli space 
approximation for the long-range dynamics. In our case the Lagrangian for the dipole interaction 
agrees with the Skyrme Lagrangian evaluated on the relativised product ansatz. 

The aim of this paper is two-fold. Specifically we want to understand the importance and nature 
of kinematic effects in two-Skyrmion dynamics. More generally we address the wider question of 
how to approximate classical soliton interactions. In this introduction we have described the 
moduli space approximation in some detail because it provides a coherent and general attempt to 
answer that question. At the end of this paper we will give a critical assessment of the moduli 
space approximation in the light of our results. 



where M 4 is Minkowski space with signature (+, — , — , — ). We denote elements of M 4 by x, with 
coordinate \i = 0, 1, 2, 3. Sometimes we also write x = (t, x). It is often useful to parametrise 
U as 



with the constraint a 2 + 7r 2 = f%. r a , a = 1,2,3, are the Pauli matrices, ir a the pion fields and 
f n is the pion decay constant. The Lagrangian density is best written in terms of the pull-back of 
the left-invariant form U'dU 



2 The Skyrme Model 



The fundamental field of Skyrme's theory is a map 

U: M 4 ^SU(2) 







(2.3) 



(2.4) 



It reads 



f 2 



■tr(L fl L») + —tr([L fl ,L v }[L»,L»}) 
■ti(R^) + ^tr([^,/a^^1)- 
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The constant e is a free parameter of the Skyrme model and it is customary to treat /„- as a free 
parameter as well. We will adjust e and f n later to fit experimental data. We fix our units by 
setting Planck's constant and the speed of light to 1. 



Note the large symmetry group of the action corresponding to Q2.5Q . It is invariant under 

Poincare group x S l O(4) chiral (2.6) 

where 

acts on U(x) via left and right multiplication with constant SU(2) matrices. The Euler-Lagrange 
equations take the form of conservation laws: the currents 

^ = ^+4^p[^[^M (2-8) 

and 

^ = R » + ^rp[ R »AR v iRA\ (2-9) 

satisfy 

djf = dpR? = 0. (2.10) 

So far we have assumed that the pions are massless. Physical pions have a small mass of about 
Rj 140MeV (we ignore the mass difference between 7r° and 7r ± ) . This can be incorporated into 
the Skyrme model by adding 

\ftm*MU + tf-2) (2.11) 

to the Lagrangian. This term explicitly breaks the chiral symmetry to 

SO(3) isospin , which acts via conjugation with a constant SU(2) element A 

U(x) ^ AU(x)A\ (2.12) 

or alternatively, in terms of the associated orthogonal matrix 

A ab = hi(T a Ar b A^) (2.13) 

and the pion field, 

n a (x) ^ A ab rt b (x). (2.14) 
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In the presence of the pion mass term the equations of motion become 

dplf = - V?" = ~ml{U - C/t). (2.15) 

It turns out that the predictions of the Skyrme model for the static properties of a single nucleon 
do not depend much on the pion mass, as long as it is small ||13|| . For the interaction of two 



Skyrme solitons, however, the pion mass has important physical consequences. Here we will 
carry out calculations mainly for the case m n = where the theory is more symmetric and the 
mathematical manipulations somewhat simpler. In certain applications we can then take into 
account the physical pion mass by general physical arguments. 

Viewed as an infinite dimensional Lagrangian system the Skyrme model's configuration space 
C is the space of maps 

U : R 3 i-> SU(2). (2.16) 
The Lagrangian L = J d 3 x C has the usual form 

L = T-V (2.17) 



where T is the kinetic energy 



f 2 

J IT 

4 i6e 2 



T = -lft I (L Q L Q )-—ti{[L Q ,L i \[L^L l }) (2.18) 



and 



V = -il tT (L i L i )--^tx([L jl L i }[L j ,L i }) (2.19) 



is the potential energy functional. To make V well defined, one needs to impose a finite energy 
requirement on elements of C. In practice one demands that the field U tends to a constant matrix, 
which one usually takes to be the identity, at spatial infinity. This allows one to compactify space 
to S 3 and regard U as an element of n 3 (S' 3 ) = Z with an associated integer degree or winding 
number. This number is absolutely conserved during time evolution and is identified with the 
baryon number B. A subset of C whose elements have a definite winding number B is denoted by 
C B . 

The imposition of the asymptotic condition also has important consequences for the symmetry 
of the theory. The chiral symmetry is broken to 5'0(3) lsospin . The euclidean group E 3 = S'0(3) spacc x 
R 3 , whose elements (Q,X) act on fields U G C via pullback: 

U(x) i-> U(g~ l (x- X)), (2.20) 

and the isopspin rotations 50(3) isospin are symmetries of the Lagrangian and map each sector Cb 
of the configuration space into itself. There is also a discrete symmetry, the combination of parity 
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operations in space and isospace. Separately these parity operations reverse the sign of the winding 
number B, but the map 

II : U{x) i — ► U\-x) (2.21) 

leaves B unchanged. Thus for each baryon number B the symmetry group of the Lagrangian 
system (C b ,L\ Cb ) is 

S = TlxE 3 xSO(3) isospin . (2.22) 



3 Dynamics of a Single Skyrmion 



We briefly review the standard approximations used in discussing single Skyrmion dynamics and 
assess their validity. 

The usual way to obtain a particular solution of the static field equations in the sector B = 1 
is to make the hedgehog ansatz 

x 

U H (x) = exp(if(r)x-r) r = \x\, x = — (3.1) 

r 

with the boundary conditions /(0) = 7r and /(oo) = to ensure that Uh has winding number 1. 
Minimising the potential energy V restricted to Skyrme fields of hedgehog type leads to a second 
order ordinary differential for /, given in ref. [pL0| , which can easily be solved numerically. In the 



following / will always stand for that solution and the corresponding field Uh will be referred to 
as the standard Skyrmion. We are particularly interested in the behaviour of / for large r, which 
is 

+ (3-2) 

where A is numerically found to be 2.16/e 2 /^ . There is strong numerical evidence that the 
standard Skyrmion is also a minimum of V among all fields in the B = 1 sector and that all 
minima in this sector can be obtained by acting on the standard Skyrmion with the symmetry 
group S. This group is 9-dimensional, but when it acts on a field of the hedgehog form spatial 
and isospin rotations are equivalent. Hence the orbit of the standard Skyrmion under S is a 6- 
dimensional manifold which we denote by Ait- This is the manifold which is used as the moduli 
space in the B = 1 sector. It is diffeomorphic to R 3 xS*0(3) and therefore its elements, which 
are called Skyrmions, are fully specified by their position and orientation. On Ai\ the potential 
energy is constant and equal to the Skyrmion's rest mass M, which can be calculated from the 
radial function /. The result is M = 73f n /e |H|. Thus the induced dynamics is determined by 
the kinetic energy, which can be found by inserting the ansatz 

A(t)U(x - X(t))A\t) (3.3) 
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into ( |2.18| ). After a lengthy calculation one finds the induced Lagrangian [|Kj 

L x = -M + Atrii f + \mX 2 , (3.4) 

where the dot denotes differentiation with respect to t. The moment of inertia A is also calculated 
in ref. |l(| (where A is denoted by A) and found to be 53.3/e 3 f n . The angular velocity a is 
usually defined via A* A = —^ol-t. In terms of the SO (3) matrix A associated to A one can also 
write A^A = ct s where Sj, i = 1,2,3 are the generators of the 5*0(3) Lie algebra satisfying 
[si, Sj] = eijkSk- In terms of a , AtrA^r = |Ack 2 . According to the equations of motion derived 
from Li a Skyrmion moves uniformly in space and rotates rigidly at constant angular velocity. 

So far the standard, approximate treatment of the B = 1 sector. How useful are the ap- 
proximations made so far? Clearly the uniformly moving but non-rotating solution is a good 
approximation to the exact solution one obtains by Lorentz-boosting a static Skyrmion, provided 
the speed is small and Lorentz contraction can be neglected. However, in this case the moduli 
space approximation is merely a complicated way of deriving something trivial. The other, rigidly 
rotating solution found in the moduli space approximation is almost certainly a poor approxima- 
tion to solutions of the full field equations. This can be seen in a number of ways. One argument, 



frequently discussed in the literature and reviewed in ref. ||11|| , is based on the observation that 
there are centrifugal effects arising from the rotation which could pull the Skyrmion apart. We 
prefer to discuss the rotating Skyrmion from another point of view which suggests that an exact 
solution exists only if ^ and that it is not of the hedgehog form at all. Our discussion is 
based on the following simple observations. Far away from the centre of the Skyrmion the field 



equations ( |2.10|) simplify to three uncoupled free scalar wave equations for the massless pion fields 

d 2 n a 

— -An a = a = 1,2, 3. (3.5) 

Furthermore the pion fields of a Skyrmion AUh{x)A^ centred at the origin behave asymptotically 
like 

^ a ~ AA— — ^- = UX^f- (3.6) 
where {e a } is the standard basis of R 3 . Defining 

p a = AnU\A- l e a (so that * a (3-7) 

we see that for sufficiently large r each of the pion fields 7r a can be approximated by a static 
solution of the scalar wave equation with a suitable (time-independent) dipole source term (we 
will discuss the charge distribution for a scalar dipole systematically in sect. 6): 

^l-A^ = -p-W 3 \x). (3.8) 
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This suggests that the asymptotic behaviour of each pion field of a rotating Skyrmion should 
be described by solutions of ( |3.8| ) with rotating dipole moments p a . More precisely, for a dipole 
moment p(t) rotating at constant angular velocity a., i.e. p(t) = p(t)xac, equation ( |3.8| ) is solved 
by 

(f)(t, x) = - — -Ux-pJt) + (cos ar + ar sin ar)x-p . (t) 

+ (ar cos ar — sin ar)x-(p(t) x a)] (3-9) 

where a — \ct\ and p» and p ± are the components of p parallel and perpendicular to ex. This 
should be compared with the asymptotic pion fields of the rigidly rotating Skyrmion 

^(t, X )^PM±. (3.10) 

Clearly the rigidly rotating solution is only adequate if r < 1/a. Furthermore we calculate, using 
standard formulae from radiation theory, that the exact solution (3.S) will radiate away its energy 
at a rate 

This suggests that a spinning Skyrmion should similarly lose energy through radiation and we can 
use the formula ( |3.11| ) to estimate the time for which the approximation of uniform rotation is 
good. To sum up, the dipole picture indicates that, for m n = 0, a rigidly and uniformly rotating 
Skyrme field can only be a good approximation to an exact solution of the Skyrme field equations 
for times less than 

^ = ^ = ^ (3.12) 
P \p\ 2 a 2 

and for distances from the centre of the Skyrmion less than 

1 

a 

Which value should we insert for the frequency a in order to obtain numerical estimates for the 
above parameters? Spin 1/2 quantum states of a single Skyrmion are models for physical baryons 
and spin 3/2 states model the A - resonance. We fix the constants f n and e as in ref. [|1(| to fit 
the nucleon and A masses 

U = 65MeV and e = 5.45, (3.14) 

so that f n is 30% smaller than its usual value of 93MeV. Thus the frequency at which a classical 
Skyrmion would have to rotate to have the proton's or neutron's angular momentum is 

aN ~ 2~/\ ~ L5 ^ r ~ 100 MeV ' 
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We then find that 

3 7 

tufa ~ j- ~ 3 x 10~ 23 s and r max ss — ss 2 fermi. (3.16) 

A classical Skyrmion spinning at the frequency oln would therefore radiate away its rotational 
energy in a very short time if m n = 0, but it could be stabilised by including the physical pion 
mass since < m n . Note, however, that in a fully quantised form of the Skyrme model the short 
classical lifetime need not necessarily be a cause for concern: a proton or neutron could only lower 
its rotational energy by emitting a particle of spin 1/2, but there is no such particle in the model. 

The value r max should be compared with the root-mean-square radius of the Skyrmion's mass 
distribution, which is r = 0.7 fermi flH| . The form of the pion fields at r > 2 fermi does not matter 



much in the calculation of the static nucleon properties in the Skyrme model, but it is crucial for 
the long-range interaction of two Skyrmions. In studying the long-range forces between spinning 
Skyrmions, and in particular in applying the Skyrme model to the nuclear two-body problem, one 
therefore needs to bear in mind that, no matter whether m n = or ^ 0, the long-range fields 
of a spinning Skyrmion are not of the hedgehog form. 

Even if we restrict ourselves to circumstances where the rigidly rotating hedgehog field is an 
adequate approximation it does not follow that a moving and spinning Skyrmion is well described 
by (|3.3j ). Rather, one should Lorentz-boost a stationary and rigidly rotating hedgehog field to 
obtain a moving and spinning Skyrmion. More precisely, suppose that the Skyrmion is moving 
with velocity u and has position X at time t. Its orientation at time t is given by an SU(2) 
matrix A, and the orientation then changes uniformly with respect to the rest frame time, which 
we denote by t'. We set t' — at the space-time point with laboratory coordinates (t,X); at a 
later laboratory time t + St, it is given by 

t' = 1 (5t-u-(x- X)) (3.17) 

where 7 = (1 — u 2 )~ 1 / 2 . Thus the field for a moving and spinning Skyrmion is 

U(t + St, x) = A{t')U H {x')A\t'), (3.18) 

where 

x'n = j(x\\ — X\\ — uSt) and x' ± = x± — X±. 

Here the suffixes || and _L denote the components parallel and perpendicular to the velocity vector 
u and A(t') is an 5'L r (2)-valued function defined by 

A(0)=A and - % -cx-r = A\t')^p-. (3.19) 

where a is a constant vector, representing the angular velocity in the Skyrmion's rest frame. 
The Skyrme field ( |3.18| ) takes into account Lorentz contraction of a moving Skyrmion and the 
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relativity of time: if the pion fields' rotation is synchronised in the Skyrmion's rest frame then 
their rotation will be retarded or advanced in the laboratory frame, depending on their distance 
from the Skyrmion's centre projected onto u. 



If the velocity u is constant the field (p. 18 ) is a good description of a spinning and moving 



Skyrmion for all laboratory times t + St. Later in this paper we will consider Skyrmions which do 
not move uniformly but whose acceleration is small. Then one needs to perform different Lorentz 
boosts for different laboratory times and the formula ( 3.18Q is only an approximation, valid for 



small increments St. This is sufficient for us, since we will only require Skyrme fields and their 
rate of change (with respect to laboratory time) at a given laboratory time t. The position vector 
X and the orientation matrix A strictly depend on t but, for ease of notation, we will not write 
this dependence explicitly. Note also that differentiation with respect to time should be carried 
out according to 

dU , s dU(t + St,x), 

In our calculations will only take into account the velocity dependence of the field ( |3.18|) to leading 
order. To this accuracy, and with the limitations discussed earlier, we will make the approximation 

t' » St - u- (x - X) (3.21) 

and describe a single moving and spinning Skyrmion by the following field: 

U{t + St, x) = A(t')U H (€)A\t') (3.22) 

where 

| = x - X - uSt + -(u-(x - X))u 

and it is implied that we use the approximation 

A(t') « A (l 2 - l -ct-r(St -u {x- X))) . (3.23) 

When evaluating the Skyrme Lagrangian on fields of this form we only retain terms which are 
at most quadratic in linear and angular velocities. Remarkably, the evaluation of the Skyrme 
Lagrangian on the field ( |3.22| ) again yields again the Lagrangian L\, but now the angular velocity 
at should be interpreted as the angular velocity in the rest frame (|3.19|). 



For later use we introduce the momenta 

J = Act and I — AAa (3.24) 

which are the conserved quantities due to the invariance of L\ under right-translations A t— > AF 
and left-translations A i— > FA by a constant 577(2) matrix F respectively. In the usual quantisation 
of a single Skyrmion J becomes the spin operator and J becomes the isospin operator. 
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4 The Product Ansatz for Moving and Spinning Skyrmions 



In this section we outline our method for calculating the asymptotic form of the Skyrme Lagrangian 
when evaluated on fields of the product form. This method applies equally to the standard product 
ansatz, which employs fields of the form (|3.3|) for the individual Skyrme fields, and to our relativised 
product ansatz, which uses the fields (|3.22|) . For definiteness, let us consider the latter. 

Assume that, at time t, one Skyrmion has orientation A e SU(2), angular velocity in its rest 
frame a , position X, and velocity u and the other Skyrmion has orientation B e SU(2), angular 
velocity in its rest frame 3, position Y and velocity v. Then we write the relativised product 
ansatz at time t + 5t close to t for the field describing this situation as 

U(t + 6t,x) = U {1) (t + St,x)U {2) (t + St,x). (4.1) 

Here U {1) is given by ( |3.22| ) and U^ 2 > analogously by 

U {2 \t + 8t,x) = B{t")U H {rj)B\t") (4.2) 

where 

t" = (1 - v 2 )-^(5t-v-(x-Y)) « St- v-(x-Y) 
T) = x — Y — vSt + -(v-(x — Y))v 

and the function B{t") is approximated by the following expression, adequate for our purposes: 

B{t") w B(l 2 - % -3-T{8t-v{x-Y))y (4.3) 

When inserting the field ( |4.1| ) into the Skyrme Lagrangian we first carry out all necessary 
differentiations and then set St = 0. For our calculations it will also be useful to introduce 

R= X -Y, R=\R\ (4.4) 

for the relative position. The computation of the Lagrangian is simplified by the observation that 

L, = U^d.U = UW(LM - Rf)uW (4.5) 

where L« = U^d^U^ and R® = U^d^U^l When inserting this expression into ( |2.5|) the 
conjugation with drops out. One still generates a large number of terms, but they can be 
further classified according to their degree in L$ or R^> . The terms which involve only either 
LW or Rj® give the sum of the Lagrangians for a single free Skyrmion. When expanding the 
remaining terms in the Lagrangian in powers of 1/R, terms which are linear in Lj-p or R^ give 
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the leading contribution. Neglecting the terms which are quadratic in both Ljp and we get 



C = -^tr(L«L^) + 

- f -ft,{RfR^) + ^tr([4 2 ), R?)[R^, Rt>]) 



+ f ^ti [L^R^ + L^R^-L^R^). (4.6) 

In terms of the mass M and the moment of inertia A of a single hedgehog Skyrmion as defined 
earlier (|3.4j) we can write the Lagrangian restricted to fields of the product form (|4.1| ) 

L 2 = -2M + l -Mu 2 + ^Aa 2 + l -Mv 2 + l -K(3 2 + L m (4.7) 



where 



L int = f -f j <?x tr (L^R^ + LfR^ - LfR^) . (4.8) 

Our aim is to calculate the leading terms in an expansion of L int in powers of 1/R. Our calculations 
are based on the idea of treating the field of one soliton as constant in the vicinity of the other. 
In order to evaluate the integral (|4.8|) we divide the region of integration into three parts 



R 3 = B p (X) U B p (Y) U V P (X, Y) (4.9) 

where B p (X) is the ball of radius p around X and 

V p (X, Y) = R 3 - (B p (X) U B p (Y)). (4.10) 

With this division we can exploit the simple asymptotic form of the hedgehog profile. Equation 
implies that, for large |£|, 

L^\t,x)r.im(t jX ):^iXd,(A(t')^(t')] (4.11) 



\6t=0 



and similarly for large \rj\ 



R?Kt,x) ~ rf\t,x) := -iXd, [B{t")^B\t")\ . (4.12) 



|<5t=0 



Since L$ and U^> differ only by a triple product of the Z/j^'s (similarly for Riy) we also have, in 
the above limits, 

L«~Z« and Rf~rf. (4.13) 

As R becomes large we will choose p so large that outside B p (X) and B p {Y) we can replace L^p 
and Ljp by as well as R^ and Rff by r^. But we also want to keep p small enough so that 

13 



ijP and may be considered constant over B P (Y) and B p (X) respectively. We achieve this by 
letting p — > oo as R —>■ oo in such a way that p/R —* 0. We then find 




(4.14) 



This formula is the starting point in most of our calculations in the following sections. 



5 The Two-Skyrmion Potential Revisited 



As a first step we consider the static situation, where u = v = and A(t') = A and B(t") = B 
are constant SU(2) matrices. Then —L int reduces to the two-Skyrmion potential V2: 

■J IT 



V 2 = ^/^tr^^+Zf^f)-^^ 



- f [tr (rPiX) J Bp{X) d*x (Z« - + tr (l?\Y) d* X (R^ - 

+ [ (FxtityprW)]. (5.1) 



/R3 

The asymptotic form of V2 was already given by Skyrme in ref. |3J and the potential energy of 
static product fields for any value of R was studied numerically in ref. [0. We will re-derive 
Skyrme's result using a different method which we find mathematically more satisfactory. In 
the course of our calculation we will also derive a number of formulae which we will need in the 
remainder of the paper. 

First we introduce the £0(3) matrix 

O ij = -tv^A^B^B^A) (5.2) 

specifying the relative orientation of the two Skyrmions. The main input of our calculation is the 
fact that both and satisfy the Euler-Lagrange equation for static fields: 

d.iP = d.R^ = 0. (5.3) 

Let us for the moment concentrate on and write it as 

LP = tAL m r m Al (5.4) 



Then, for each the components of a smooth vector field on R 3 and ( |5.3|) 

implies that this vector field is divergenceless. Thus, for each m we can find a vector field with 
components Z^m such that 

Lira = tijkdjZkm- (5-5) 
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Explicitly we find 



rr i v \ ( ■ 2 j-/ \ 1 sin 4 /(r)\ _ (l ... N 1 ,sin 2 /(r)\ 
Z fcm (a; + A) = sin f{r) — - x k x m + -rf (r) + —^f (r) e kmn x n . 



(5.6) 

Hence 

j d?xLf\x) = iAr m A ] [ dS e ijk Z km (x + X)xj 

JBJX) JdBJQ) 



,{X) JdB p (0) 

iAr m A^ I dS {{8 im - XlXm )(\rf\r) + J-/' (r) — /(r) 

iAnA+y (Vf (p) + ^P/'(P) sm 2 /(p)) 
8?rA ^^t. (5.7) 



3 

In the last line we have used that p is so large that we can use the asymptotic expression (|3.2| ) for 
/ and have omitted non-leading terms. Combining this with 

/ v d?xlf\x) = -i\Ar m A^ [ v d 3 xd t d m . 1 
Jb p (X) Jb p (X) \x — X\ 

lr 1 

= -i\Ar m A^-5 im d 3 xA— 
3 JBJO) \x\ 



we have 



Since 



(o) 

-—Ar t A j (5.* 



[ d 3 x - lP) ~ -47rAzA^A t . (5.9) 



rf } (X) = -zA — ^ Rn Br n B^ (5.10) 



we therefore find 

tr (rf»(X) A (!<" - i<»)) ~ -8, / r0 ^f 0B (5.11) 

An analogous calculation for the integral over B p (Y) yields also 

tr (l?\Y) d"x (R? - rf')) ~ ^x' tT°Z^°t . (5.12) 

To calculate the remaining integral we integrate by parts and use again A(|;r| _1 ) = — 4n5^(x): 

{ d 3 xtr(A 2) lp) = 2\ 2 O lk ( d 3 xdid ri — 1 -—d l d k - l — 1 —— 
7rs \ 1 1 / J R s J \x — X\ \x — Y\ 

= 8ttA 2 — . (5.13) 
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Combining terms and introducing the constant 

k = 2v\ 2 fl (5.14) 

we obtain the result: 

v 2 „ _ 2K t r0 -3(A. A) = i 

it' 3 it 

Here V is the gradient operator with respect to R. In the following we will use the symbols d{ and 
V for differentiation with respect to either x or R. Which of the two possibilities applies will be 
clear from the context. Note that one would get — V2 from integrating the asymptotic expression 
tr f/| 1 V| 2 ^ over the whole of R 3 . The effect of replacing the singular part of and rf^ by the 
true, smooth solutions inside the balls is to add 2x V2. 

One can usefully write Vj> by expressing the rotation matrix O in terms of its axis n and the 
angle of rotation about that axis tp, i.e. 

Oij = cosipSy + (1 — cosip)nirij + sin ^e^n^. (5.16) 

and finds 

V 2 „ _ 2/c (i _ cos^) 1 " 3 ^'^ 2 . (5.17) 

Although the details of our calculations may seem quite complicated, the essential steps are 
simple. After suitably dividing the region of integration we were able to express the integral for 
involving U(x) at all points in R 3 in the simple form ( 5.15Q which requires only information 



about the asymptotic nature of the fields. This remarkable result rests essentially on the fact 
that the two B = 1 Skyrmions that enter the product ansatz individually satisfy the equations of 
motion, as expressed in (|5.3|). 



There is another, heuristic way of deriving ( 5.15Q . It is based on the observation, explained in 



sect. 3, that, from afar, a Skyrmion looks like a triplet of orthogonal scalar dipoles. Thus it is 
plausible that the long-range static potential between two Skyrmions can be understood in terms 
of scalar dipole interactions. Specifically one considers two triplets {p a } and {q a }, a = 1,2,3 of 
scalar dipoles. There is a scalar dipole interaction between p a and q b only if a = b, given by 

i-(p o .V)(q a -V)i (5.18) 

To obtain the total potential we simply sum the interactions of the pairs. We also need the formula 
for the dipole moments associated to a Skyrmion in general orientation ( |3.7| ). Let A and B be the 
SO (3) matrices associated to the SU{2) matrices A and B as in ( |2.13| ). Then we have from ( |3.7| ) 

p a = Anf^XA^ea and q a = Anf n \B- l e a (5.19) 
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so that p a = Oq a . Thus the potential for the interaction of two triplets of scalar dipoles is 

V 2 = 2k(V-OV)4 (5-20) 
which agrees with the asymptotic form of Vi- 

6 Skyrmions as Triplets of Scalar Dipoles 

In this section we take the dipole model further: we calculate the velocity-dependent forces between 
two triplets of scalar dipoles. 

How do moving scalar dipoles interact? To answer this question we need to know the coupling 
between a massless scalar field and a general source and to write down the charge distribution for 
a moving and spinning scalar dipole. We fix the coupling of a scalar field to a general charge 
distribution p by writing the field Lagrangian density as 

Aeld = 7^>9> + P <t>. (6.1) 

For given p this leads to the wave equation 

d 2 (j) 



dt 2 

which is solved by 



P (6-2) 



l r * ,p(t — \x — x'\,x') 

4>(t,x) = — / d 3 x- — r r—*-- 6 - 3 

v ' Air J \x-x'\ 

Turning to the charge distribution p we notice that it must be a Lorentz scalar to make equation 
( |6.2| ) consistent. As a first step we give the charge distribution of a moving point particle. In its 
rest frame is is characterised by a coupling strength g and its position X and is given by 

g5 i3 \x-X). (6.4) 

To find its charge distribution in a general frame we write down a charge distribution which 
is manifestly a Lorentz scalar and which reduces to (|6.4j ) in the particle's rest frame. For a 
particle with world line X{r) = (T(r), X(r)) and four-velocity u = (u°, u) = dX/dr the required 
distribution is 



/W (x)=gj dr ^~^5 {4 \x - X{r)). (6.5) 
Carrying out the integration over the world line parameter r and using 

(^Pj = Vl - w 2 (6.6) 
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we find 



Ppoint (t, x) = gy/T=tf5®(x - X(t)) (6.7) 

where, in an abuse of notation, we have written X(t) for X(T~ 1 (t)). 

Our model for a scalar dipole involves a world line X(t) and the dipole moment p(r) in the 
dipole's rest frame. In that frame, which is characterised by u = dX/dr = (1,0), we can think 
of the dipole as consisting of two equal and opposite scalar point charges with rest charges ±g , 
situated at X + \d and X — see fig. 1. 




Fig. 1 



Dipole charge distribution in the dipole's rest frame 



Still in the rest frame we impose that |d| remains constant during time evolution and introduce 
the angular velocity ol 

dd 

— = dxa. (6.8) 
dr 
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Then the charge distribution depicted in fig. 1 is described by 



9 \ 1 " [at) 6(3){X ~ {X + \ d)) ~ 9 \ l - 6(3){X - {X ~ \ d)) - (6 ' 9) 
In the limit where \d\ becomes small and g large in such a way that the dipole moment 

p = gd (6.10) 

remains constant this simplifies to 

-p-V5 i3) (x- X). (6.11) 

Again we write this charge distribution in a general frame by giving a Lorentz scalar distribution 
which reduces to ( |6.11| ) in the rest frame u = (1, 0): 



p dipol e(x) = - J dr ^~^P,(r)d^ (4 \x - X(r)). 



(6.12) 



Here the four-vector P reduces to the dipole moment p in the dipole's rest frame, so in particular 
u^PV 1 = 0. This condition determines P in terms of p: 



P° 
P 



t(«-p) 



P + 



7 



7 + 1 



u-p)u 



(6.13) 



where 7 = (1 — u 2 )- 1 ' 2 as before. Making these replacements in (|6.12|) and carrying out the 
integration over the world line parameter r we arrive at the expression 



/ \ r d , \ 1 r-, ( 7 1 
Pdi P oic {t,x) = -[—(u-p) + -p-V+ —— - - 
at 7 V 7 + 1 7 



{u-p)(u-\7)}5 3 (x- X(t)). 



(6.14) 



Here p = p(T 1 (t)) should be thought of as a function of time. Assuming du/dt ~ and retaining 
only terms which are at most quadratic in the velocities u and cx this reduces to the distribution 

1 _ 1 



Pl (t, x) = -[(u-pxa) + (1 - ^-u 2 )p-V - -(u-p)(u-V)]5 3 (x - X(t)). 



(6.15) 



We can now calculate the potential due to p\ according to (|6.3|) . Since p t varies slowly with 
time for slowly moving sources we can make some approximations: 



(f>i(t,x) 



1 

47T 



, 3 ^,pi(t,x') d 



d x 



\x — X 



dt 



J d 3 x' Pl (t,x') + J d 3 x' Pl (t,x')\x' -x\\ 



(6.16) 



We should be careful to understand the range of validity of our approximation: for constant p 
but time-varying X we can neglect higher order terms in the above Taylor expansion if u is small 
relative to the speed of light and the acceleration and higher derivatives are small. For fixed X 
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and time dependent p, however, higher terms in the Taylor series are not necessarily small even 
if p rotates uniformly In that case the dipole field is exactly given by ( |3.9| ) and ( |6.16| ) is an 
expansion in powers of |o;||aj — X\ of the exact solution. Retaining only terms that contain at 
most two time derivatives we obtain 

0i(t, x) « i- {{u-cxxp) - (1 - l -u 2 )p- V + I( w . p )( w .V))^4x| ~ \lp P ' V|£C " X| ) • 

(6.17) 

Now consider a second scalar dipole with world line Y{p) = (S(a),Y(a)) and with dipole 
moment q in its rest frame. Denote its three- velocity dY/dt by v and its angular velocity in its 
rest frame by (3, i.e. dq/da = qxf3. Its charge distribution, in the same approximation as for the 
first dipole, is 

p 2 (t,x) = -[(v-qx/3) + (1 - l -v 2 )(q-V) - ±(vq)(vV)]6 3 (x - Y(t)). (6.18) 

The part of the Lagrangian which describes the interaction of the second dipole with the field 
produced by the first is therefore 



L int = J d 3 x'p 2 (t,x')(j) 1 (t,x') 



~ —[(V-p)(v-/3xq)-(q-V)(u-ocxp) 

~ {l-\u 2 - l -v 2 ){p.V){q.V) 

+ \{u-p){u-V){q-V) + \{v-q){v-V){p-V)]^ 

- i^ v 4^ K (6 - 19) 



As in the previous section R = X — Y , R = \R\. In the last line V differentiates with respect to 
R and d/dt, when operating on R, acts on X only. It is useful to rewrite the last term 

(q-V)- 2 ( P .V)R = - [( q -V) Jt (p-V)R) 

~(q-V)(p-V)R- (q-V)(p-V)(u-V)R 
+(q-V)(p-V)(v-V)R+(q-V)(p-V)(u-V)(v-V)R. 

(6.20) 

The time derivative outside the big brackets in the first line acts on all the dynamical variables 
P,q,X,Y so 

I (( 9 -V)|(fV)«) (6.21) 
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is a total time derivative which we may discard since we are working with a Lagrangian. 

Finally we translate our results into "Skyrme language" , using the notation introduced in sect. 
5. Discarding the total time derivative in ( |6.20| ) and summing over three pairs of dipoles p a and 
q a we find, using the formulae ( |5.19|) , 



L int = k[2(V-C?(ux/3)) - 2(nxa-CV) 

+ (u 2 + v 2 )(y-OV) + («-V)(u-OV) + (u-V)(V-C?«)]-^ 

+ k[(Vxq-0(Vxj9)) 

+ («-V)(Vxa-OV)-(u-V)(V-0(Vx/3)) 

- (V-OV)(u-V)(v-V)}R 

- 2/c(V-0V)— . (6-22) 

R 

Note that the terms in this Lagrangian can be classified into three types: spin-spin couplings 
which couple the angular velocities a and f3 to each other with strength proportional to 1/R, 
spin-orbit couplings which couple a and (3 to u and v with strength proportional to 1/R 2 and 
orbit-orbit couplings which are quadratic in u and v and which are of order 1/R 3 . Before we 
discuss our result further we derive it from the relativised product ansatz. 



7 Two-Skyrmion Dynamics according to the Relativised Product 

Ansatz 



In this section we evaluate the interaction Lagrangian L int on the time-dependent relativised prod- 
uct ansatz up to and including terms which are quadratic in the velocities u,v, a and (3. For 
the sake of the impatient reader, who may wish to proceed immediately to the comments at the 
end of this section, we announce already at this point that the asymptotic form of our final result 
agrees, up to a total time derivative, with the interaction Lagrangian L int calculated in the previous 
section. 

We begin with the following kinetic terms of the interaction Lagrangian: 

T int = | / d 3 *tr (L«4 2) + W - • (7-1) 

For time dependent fields, L int actually contains further terms with explicit time derivatives in the 
double commutator [L^\ [Lq,L^]\ contained in uf 3 (similarly for Rf^) but for computational 
purposes it is best initially to concentrate on T int . Since each term in T int explicitly contains 
two time derivatives we can omit the relativistic corrections in the product ansatz (|4.1|) when 
calculating ([7.1|) since we only retain terms quadratic in the velocities anyway. Thus we replace 
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t' and t" by 5t, £ by x — X — u6t and r] by x — Y — -u^t in ( |4.1| ). To carry out the integration 
in we employ the same technique as for the static terms, so we divide R 3 as in Q4.9Q , extend 
the integration for the asymptotic expressions of the fields from V P (X, Y) to the whole of R 3 , 
subtract the integrals of the singular parts over the balls and add the integrals of the true fields 
over the balls. This time we will omit most of the details, but it is instructive to write down the 
integrals of the asymptotic expressions explicitly. These now read 

l { o\t, x) = iXA UijkCtiTjdk ^ x + UiTjdidj - — — — j A 1 (7.2) 



and 



so that 



rf\t, x) = -i\B (eijkPiTjdk- — l — - + v^didj- — l —- ) B ] (7.3) 



\x — r \x 



f d 3 xtr (/o 1} ^ 2) ) = f^\ 2 [e ijk e lmn O jm aiPi f d 3 xd k - l -— d n - — - 

jR^ J \X — y\. \X — 

f 1 1 

+ e ijk aiViO jm / d 3 x d k - — did m - 

J \x — X\ \x — 



Y\ 
~Y\ 



f 11 

+ ClmnPlUiOjrn d 3 X didj- —r 9 n - 

J \x — X\ \x — 

f 1 1 

+ UiViO jm d 3 xdidj- — -did m - — ]. (7.4) 

J \x — X\ \x — Y\ 



We use the first of the above integrals to illustrate our method of integration: integrating by parts, 
shifting the origin to Y we get 

d 3 x d k - 3—— d n - — ^ . = — [ d 3 x -r^-r d n d, 



I | i35 — JT | J \X\ \X — K\ 

_ f 3 I d d 1 

; J \x\dR n dR k \x- R\ { } 

Next we would like to take the differentiation outside the integral. To justify this interchange of 
limits we consider the above integrals over a large , but finite ball -8^(0) with radius K > R, 
centered at the origin. Then 

d 3 x- — - — ~ = AttK -2tiR. (7.6) 

B K {0) \X — it | \X\ 

For the finite region of integration, the exchange of limits is allowed, so 

d 3 x — - — - — r = - — - — Ukk-2iiR) 

b k (o) \x\ oR n oR k \x — R\ oR n oR k 



-2nd k d n R. (7.7) 
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We see that, as long as K > R, the integral on the left is actually independent of K, we can take 



the limit K — > oo and recover the original integral. Evaluating the other integrals in ( |7.4j ) in a 
similar fashion we obtain 



^ f d 3 xtr(l { Q } 4 2) ) = K[(Vxa-0(Vx(3))R 

+ («-V)(Vxa-0V)i2 



(u-V)(V-0(Vx/3))i2 

(V-OV)(u-V)(v-V)R}. (7i 



the static case, this amounts to adding 



How is this result affected when we take into account the true form of L and Rq ? As in 



hr(ri?\x) f x/tiW-W)) (7-9) 



and 



|tr^(r)/^ )C / 3 x(4 2) -rS 2) )). (7.10) 

The net effect of this is to change the signs of some but not all of the terms generated when the 
differentiation in f[7.8|) is carried out. The terms in the integrand of (|7.9|) which are proportional 
to a and the terms in the integrand of Q7.10 ) which are proportional to f3 lead to spin-spin terms 



which fall off like 1/R and to spin-orbit terms which fall off like 1/R . Such terms are non-leading 
and we neglect them. Thus for our purposes = —mLf^ and R^ = —ViRf\ Hence we can 
again exploit the fact that the vector fields L and R are divergenceless, but we do not give 
the details here. The result is 

T int ~ K[(Vxa-0(Vx(3))R 

+ (v-V)(Vxa-0V)i?-2(Vxa-CV)4 

R 

- (u.V)(V-0(Vx/3))i*+2(u.0(Vx/3))- 

R 

- (V-OV)(u-V)(v-V) J R+2(w-OV)(v-V)4+2(V-Ov)(w-V)4]- 

R R 

(7.11) 



Next we evaluate the potential energy V2 ( |5.1| ) on the relativised product ansatz (fOF), again 
setting 5t = after the calculation. The velocity dependence of that field will lead to velocity- 
dependent modifications of the potential calculated in sect. 5. As a first step we restrict attention 
to constant matrix functions A(t') and B(t"), i.e. vanishing angular velocities. We then introduce 
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a "rest frame current" with components via 

-C 



L W = t7 (Dt^(i) (712 ) 



and find 



L« = Lg J + -u^Lfh (7.13) 



Further we define the current with components 
which satisfies 



|4f - 0. (7,5) 



d_ 

dr], 

and analogously to L9^ we also have 



(7.14) 



Carefully taking into account the explicit time derivatives in the double commutator 
[L ( q\ [L ( o\lP]] contained in jjf* one checks that 

Defining 



i2g) = U^i-UW (7.17) 



R? ] =R^ + \v l {v 3 m. (7.18) 



2 tv J 

Then, when integrating L^(t, x) over B p (X), we express in terms of according to ( |7.13| ) 
and change integration variables via d 3 x = (1 — ^u 2 )d 3 ^. Thus, to order u 2 , 

[ #xL<p(t,x) = ((1 - \u 2 )5 m + l -u iUm ) f d 3 afl(t,t) 

= -((1 - ^u 2 )5 im + ^UiU m )^^AT m A^ (7.19) 

where we used ( |5.7| ) in the last line. The integral of Rf\t,x) over B p (Y) can be dealt with 
in an analogous fashion. To calculate the remaining integrals in Vi we need expressions for the 
asymptotic form of and 



lf\t, x) = -iXAr m A^[[{l - -u 2 )^ - l -{u-V)u m d^j * + ^(u-V)(u-V)did m \x - X\] 



(7.20) 
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and 

rf\t, x) = iXBr m B^[[{l - ^v 2 )^ - ~{v-V)v m a^ ^ + ^{v-V){v-V)did m \x - Y\}. 

(7.21) 

Using these formulae it is not difficult to reduce the integrations of l!p , r| and tr(/| 1 V^) to 
integrals of the same form as those appearing in T int earlier in his section. The resulting asymptotic 
expression for V2 is 

V 2 ~ 2*[(1 - iti a - ^ 2 )(V-0V) + i(u-V)(«.OV) + I(«.V)(V.O«)]i (7.22) 

We still have to allow for A(t') and B(t") to vary with t' and t" respectively, which leads to 
both being a function of x. This produces an additional term in the formula (f7.20|) for if' . 

iXAdi (u-(x — -X") (axV-r) ^l )^ ^' 2 ^ 



and a similar extra term in the formula (|7.21|) for 



,(2). 



iXBdi\v(x-Y)(BxV-T)-. — Kttt I B ] . (7.24) 
V kc - Y I 



Inserting these extra terms into (5.1) one thus expects to obtain extra spin-orbit terms of order 
1/R 2 . Further spin-orbit terms are generated if one considers the effect of the relativistic time 
transformation on and R\ , but these are of order 1/R 3 and non-leading. Miraculously, 

however, the sum of the extra spin-orbit terms generated when the above corrections to and 

(21 

r\ are taken into account vanishes. Since all the integrals required to prove this result have 
appeared earlier in this section we will not go through the details here. Physically this result 
means that the relativity of time does not affect the long-range interaction between two Skyrmions 
to the accuracy considered here: we could have replaced t' and t" by St straight away. 

Before we combine T int and V2 into our final expression for L int we exploit the fact that we can 
discard total derivatives in a Lagrangian. Noting that, for any two vectors r and s 

(r-Os) = (ctxr-Os) + (r-0((3xs)) (7.25) 

and neglecting terms proportional to u (which are of order 1 / R A ) we find that the term 

[(u-0(Vx(3)) + ( w .0V)KV)]i, (7.26) 

which occurs in T int , can be replaced by the term 

[(axw-OV) + (w-OV)(w-V)]- (7.27) 

R 
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because the two terms differ only by the total derivative 4((ti-OV)i2 1 ). Similarly we can replace 



by 



-(Vxa-Ov) + {V-Ov){u-V)}- (7.28) 

R 



[(V-0(«x/3)) + (V-0v)(v-V)]4 (7-29) 

R 



since the two terms differ by 4 ((V-Ov)i? x ). 

With these replacements we have finally shown that the asymptotic form of L iat is 

L int ~ L int . (7.30) 

We already classified the terms in L int at the end of the previous section, but the result ( |7.30| ) 
deserves some further comments. First note that the Lagrangian 

L 2 = -2M + X -Mu 2 + hvlv 2 + l -k(x 2 + + L int (7.31) 

is not the asymptotic Lagrangian which one would obtain in the usual moduli space approximation: 
we have evaluated the Skyrme Lagrangian L ( p. 17] ) on a set of velocity-dependent fields which, at 



a given time, differ from fields in the moduli space M 12 by terms of order u 2 and v 2 . As a result 
the velocity-dependent terms in L 2 are not simply the restriction of the kinetic energy T to the 
submanifold M 12 of C but also include contributions from the potential energy evaluated on the 
relativised, and hence velocity-dependent product ansatz. 

Our second comment concerns the validity of L 2 as an approximate Lagrangian for two- 
Skyrmion dynamics. This is restricted in two ways. On the one hand it is well-known that 
the product ansatz is only a good description of two well-separated Skyrmions. To qualify "well- 



separated" we note that the numerical work in ref. |15[ indicates that the product ansatz is already 
a good approximation for Skyrmion separations of R > 1/2 « 1.5 fermi. On the other hand the 
allowed range of R is restricted from above for spinning Skyrmions: for given individual angular 
velocities ex and (3 the field of one Skyrmion at the position of the other is only well-described by 
the rigidly rotating field implicit in the product ansatz if 

\a\R<l and \/3\R<l. (7.32) 

This assumption was made explicitly in our dipole calculation. When the individual Skyrmions are 
spinning at the nucleon frequency a at (|3.15|) the condition ( 7.32 ) becomes R < 2 fermi. In terms 



of the standard classification of nuclear forces according to their range |16| this restriction means 
that the product ansatz, even in its relativised form, can only be used to calculate the Skyrme 
model's predictions for intermediate-range nuclear forces. This problem is routinely ignored in the 
literature; we will come back to it in sect. 9. 
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Recall, finally, that effects due to the relativity of time turned out to be irrelevant, at least 
to the accuracy considered here. This enabled us to replace the rest frame times t' and t" by the 
laboratory time St. Thus we can, at a given time t, define the relative orientation O = A~ l B, 
which in turn is essential for the notion of attractive channel fields, to be considered in the next 
section. 



8 Dynamics in the Most Attractive Channel 

The truncated Lagrangian system with Lagrangian L 2 is still rather complicated. To understand 
it better we restrict attention to a particular set of fields and truncate the Lagrangian system 
further. 



It is clear from (|5.17 ) that the static two-Skyrmion potential is most attractive if one Skyrmion 



is rotated by 180° relative to the other about an axis orthogonal to the line joining the two 
Skyrmions. Such fields are said to be in the most attractive channel and they can all be obtained 
by acting with the symmetry group S (|2.22|) on a family of standard fields, parametrised by the 
separation parameter R. In terms of a cartesian basis {ei, e2, 63} of physical space such a standard 
field consists of one standard hedgehog at \Re 3 and another hedgehog, rotated by 180° about e 2 , 
at — |i?e 3 : 

U R {x) = U H (x - ^Re 3 )r 2 U H (x + ^Re 3 )r 2 . (8.1) 

We can write the action of the group E 3 x SO(3) isospin G S in terms of an SU (2) matrix G and its 
associated 5*0(3) matrix Q for spatial rotations, a translation vector S, and an 577(2) matrix C 
for isorotations. Acting with (G, S, C) on Ur{x) we obtain 

C&U H {x -S- X -RQe 3 )Gr 2 G^XJ H {x -S + ^Rge 3 )Gr 2 Cl (8.2) 

Comparing this expression with our previous parametrisation of product fields we find the following 
relationship between the set {A, B, X, Y} and the set {G, S, C}: 

A = CG^ B = -iCt 2 G ] 
X = S + ^RQe 3 Y = S- ^RQe 3 . (8.3) 

Thus in particular R = Qe 3 . When G, S and C vary with time we also need notation for the 
associated velocities. For spatial rotations we define the space-fixed (or right-invariant) angular 
velocity 7 via 

- l - 1 . T = G& (8.4) 
as well as the body-fixed (or left-invariant) angular velocity uj 

- 1 -u-t = G ] G (8.5) 
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which are related via Quj = 7. Introducing the left-invariant angular velocity fi in isospace 

-Ul-T = C^C (8.6) 

we find the following relationships: 

a = g(n-u>) /3 = g(v 2 n-uj) (8.7) 

where T>i is the SO (3) matrix for a rotation by 180° about the z-th axis. For the linear velocities 
one checks that 

u = S + ^(RR + jxR) v = S -~(RR + jxR). (8.8) 

Using these relationships we shall now express the restriction of L 2 to relativised product fields 
in the most attractive channel in terms of u>, f2, R, S and R. This requires lengthy calculations 
which we will not describe in detail. Note, however, that it is convenient to use 

Oij = -5 i:i + 2mnj, (8.9) 

where n = Qe 2 , to simplify the terms in L 2 . One finds for example 

(«-V)(«.V)(V-OV)i2 = -4 (n ' n X^ n) . (810) 

R 6 

After simplifying L 2 in this fashion it is straightforward to replace the velocities a, (3, u and v as 
in Q and 

Before we state the result we note some features which one can derive without any calculations. 
Combining the inversion map n (|2~21| ) with rotations by 180° in both space and isospace we obtain 
9 isospace-space reflections 

U ai : U{x) ^ T a U\YiiX)T a (8.11) 

where Ilj is the reflection in the plane orthogonal to ej. Writing these maps more explicitly in 
terms of the pion fields we have for example 

Ili3 : 7r(cc) h-> (-tti(xi, x 2 , -x 3 ), 7r 2 (xi, x 2 , -x 3 ), ir 3 (xi, x 2 , -x 3 ))- (8.12) 

The maps Il a j are symmetries of the Skyrme Lagrangian and, since they leave the set of product 
fields invariant, they are also symmetries of the truncated Lagrangian system with Lagrangian L 2 . 
In order to understand the action of the maps Il a j on the attractive channel fields it is sufficient 



to consider a standard field Ur(x) (|8.1|). Such a field is only strictly invariant under the action 
of n 13 but it also invariant under the action of LTn and U 22 modulo the commutator of two well- 
separated hedgehog fields. In our calculations of the asymptotic forces such commutators are 
ignored. This can be seen from the fact that L 2 is symmetric under the exchange of the individual 
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Skyrmions position, orientation and linear and angular velocities, although the product ansatz 
(|4. 1|) is not invariant under these operations. It was shown by Verbaarschot et al. JTJJ that the 
true attractive channel fields in standard orientation and position, obtained by minimising the 
static potential energy for a fixed value of R, are in fact invariant under the symmetries n 13 , IIn 
and IT 22 . This suggests that one should only trust the attractive channel fields of the product form 
in the asymptotic limit, where they also have those symmetries. 

The invariance (modulo commutator) of Ur(x) under II13, IIn and H22 implies that, for each 
R, these maps act on the orbit of Ur{x) under E 3 x SO(3) isospin . Explicitly we find, for (a,i) = 
(1,3), (1,1), (2, 2): 

U ai :S^U t S G^nGn C^r a Gr a . (8.13) 

For the angular velocities this implies 

U ai : ui \-> V t u Ct h+ V a n. (8.14) 

The restriction of L 2 to the orbit of the one-parameter family of fields {Ur(x)} under £3 x 
SO(3) lsospln , which we denote by L att , must be invariant under ( |8.13| ) and (|8.14|) and also under 
the left action of E 3 x SO(3) isospin on itself. This implies that the most general form of L att in the 
centre-of-mass frame (S = 0) is 

/xi? 2 + a 2 uf + d 2 cu 2 + c 2 ul + A 2 ttj + B 2 Vi\ + C 2 Vl\ + 2eu 2 tt 2 + ^ (8.15) 

where a, c, d, e, A, B, C and /x are functions of R only. For later use it is best to "complete the 
square" and to rewrite this expression slightly: 

fxR 2 + a 2 ul + b 2 u 2 2 + c 2 ul + A 2 nj + B 2 (tl 2 + wcu 2 ) 2 + C 2 nj + (8.16) 

R 6 

where b 2 = d 2 — e 2 /B 2 and w = e/B 2 . 

Carrying out the explicit calculation of L att we find 

,M 2 15k. 2 ,M„ 2 a 7k. 2 2 

+ ( T ^ 2 + A-— V 2 + ( T i? 2 + A-— V 2 2 + A^ 3 2 

+ (A + + (A - + Anj - 2(A - |> 2 fi 2 + ^. (8.17) 

Note that this Lagrangian is not Galilei invariant. Although the momentum Ps conjugate to S is 
conserved one finds that in general S 7^ 0. We can, however, set Ps = and it then follows that 
S = 0. L att then takes the form ( |3.15|) and we read off the coefficient functions /x, a 2 ,c 2 ,A 2 ,S 2 
and C 2 directly; after bringing it into the form ( B.16j ) we also find 

b\R) = fR 2 -^ and w(R) = -(1 - ^) (8.18) 
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with corrections of order 1/R 2 . 

We remark that the kinetic part of L att can formally be interpreted in terms of a Riemannian 
structure on the space of attractive channel fields. Such a geometric interpretation plays an 
important role in the moduli space approximation to the dynamics of Bogomol'nyi-type solitons. 
In that context the geometric viewpoint is emphasised because the kinetic energy of the field 
Lagrangian can itself be thought of in terms of a Riemannian structure on the configuration space 
C. It is then consistent to interpret the kinetic energy of the truncated Lagrangian in terms of an 
induced Riemannian structure. Recall, however, that in our case the kinetic part of the Lagrangian 
L 2 is not just the restriction of the kinetic energy ( [2.18|) to some submanifold of the configuration 
space C. Thus the geometric viewpoint, though possibly interesting, seems less well founded in 
our case and we will not emphasise it in the following. 

The asymptotic form of the functions A 2 , B 2 , C 2 

A 2 ~B 2 ~C 2 ~ A (8.19) 

can easily be interpreted in terms of the Skyrmions' moments of inertia. With our normalisation 
of the kinetic energy the formula (|8.19|) implies that the asymptotic moment of inertia in isospace 



is 2A about any axis, which is what one expects for the simultaneous isorotation of two Skyrmions 
with moment of inertia A each. Similarly the asymptotic behaviour 

a 2 ~b 2 ~ ^-R 2 + A, c 2 ~ A (8.20) 

can be understood in terms of the moments of inertia for the rigid rotation of two Skyrmions 
a distance R apart: 2A is the moment of inertia about the line joining the Skyrmions and, by 
Steiner's theorem, (2M)(-|) 2 + 2A is the moment of inertia about either of the two perpendicular 
directions through the midpoint of the line joining the Skyrmions and at right angles to it. 

The corrections to the asymptotic expressions ( |8.19| ) and ( |8.20| ) reflect the velocity-dependent 
forces between moving and spinning Skyrmions. To study their effects on two-Skyrmion dynamics 
we will now briefly investigate the equations of motion corresponding to L att , using the Hamiltonian 
formalism. Thus we introduce body-fixed angular momenta in space and isospace 

m = — — M = — — 8.21 

OOJ Oil 

which are conjugate to to and O respectively; in components 

m x = 2o 2 lu 1 , m 2 = 2(b 2 + B 2 w 2 )uo 2 + 2B 2 VL 2) m 3 = 2c 2 u; 3 

Mi = 2A 2 VL ll M 2 = 2B 2 (Q 2 + wu 2 ), M 3 = 2C 2 Q 3 . (8.22) 

Finally defining Pr = dL^/dR we arrive at the Hamiltonian in the centre-of-mass frame 

rr Pr i rn\ (m 2 -wM 2 ) 2 m 2 M 2 M| M| _ 4k 
att 4/i + 4a 2 + Ab 2 + 4c 2 + AA 2 + AB 2 + 4C 2 i? 3 ' 1 } 
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The conserved quanitites are the Hamiltonian itself, the squared spatial angular momentum 

m 2 = m^ + m2 + ml (8.24) 

and the squared isospin 

M 2 = Ml + Mf + Ml (8.25) 

as well as the space-fixed angular momenta in space and isospace which will, however, not concern 
us here. Using the basic Poisson brackets 

{m*, rrij} = -e ijk m k , {M { , Mj} = -e ijk M k , {m,, Mj} = (8.26) 

it is not difficult to write down the Hamiltonian equations of motion F = {H, F} for any function 
F of the coordinates and momenta. Exploiting the conservation laws one can also study certain 
trajectories qualitatively without integrating the equations of motion. 

One can, for example, satisfy the equations of motion by setting S = 0, m = (m°, 0, 0) and 
M = (M®, 0, 0), where m\ and M° are two positive constants, and then determine R(t) from the 
energy conservation law 

( M_2k - K) 2 , W) 2 _ An 

l 4 R 3> MR 2 + 4A- 30k/ R 4A + 4k/R R 3 ' l ' ' 

Using the numerical values for M and A given in sect. 3 we can rewrite this equation in terms of 
the dimensionless quantity R = f n R, and simplify it, assuming R > 1: 

4£ (M?) 2 • 1 / (m?) 2 (Mj 3 ) 2 1.8\ 

-M~ — = R + ^{-R^ R-'W)- (8 - 28) 

One understands the time evolution of R(t) qualitatively by looking at the "effective potential" 
in round brackets. Remarkably this potential includes an isospin-dependent attractive Coulomb 
potential which could lead to a bound motion. For such a motion one finds that R > 2.3/M°. 
However, the angular velocity Q± = M°/2A leads, via (|8.7[), to individual angular velocities a\ 
and (3\ which are also of the order M{'/2A (one checks that uj\ is negligible) and which, by ( |7.32|) , 



restrict the validity of L att to R < 0.7/M°. Thus the bound motion predicted by equation (|8.28|) 
is always outside the range of validity of L att . 

Further studies of trajectories indicate that the result of our example calculation is generic: 
the velocity-dependent forces encoded in coefficient functions of L att give only small corrections to 
the equations of motion obtained from the following, simpler Lagrangian in which the velocity- 
dependent forces are ignored. 

L M = 2M5 2 + ^ 2 + (^ 2 + AV 2 + (^ 2 + AV 2 + Ac 2 

+ A(fi 2 + fi 2 + fi 2 )-2Acu 2 fi 2 + ^. (8.29) 
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Trajectories calculated with L att are only qualitatively different from those calculated with L stat if 
the momenta m and M are taken to be so large that L att is no longer a good approximation for 
two-Skyrmion dynamics. 

The calculations in this section and the previous one were based on moduli space ideas, with 
some relativistic corrections. The dipole picture suggests that this approach is fundamentally lim- 
ited by the condition ( |7.32|) and consequently only capable of calculating small velocity-dependent 
effects. However, the dipole model also provides an alternative method for studying the long-range 
forces between spinning Skyrmions which is in principle valid for any value of the Skyrmions' an- 
gular velocities and any (sufficiently large) value of R. This point of view is elaborated in the 
following section. 



9 Long-range Forces between Spinning Skyrmions and Nuclear 

Forces 



The semiclassical angular speed ( |3.15|) at which a Skyrmion has to rotate to have a nucleon's 
spin is so large that the product ansatz (even in the relativised form) is not a reliable basis for 
calculating the Skyrme model's prediction for the long-range nuclear force. We will therefore give 
a heuristic derivation of the long-range force between stationary but rapidly spinning Skyrmions 
based on the dipole model. 

Consider two triplets of well-separated and rapidly spinning scalar dipoles. The dipole moments 
orthogonal to the direction of the angular velocity produce a rapidly oscillating field but the dipole 
moment in the direction of the angular velocity produces the simple field of a static dipole, cf. 
(|3.7|) . When two spinning triplets of dipoles interact the oscillating fields produce no net force 
but the projections of the dipole moments onto the direction of the angular velocities interact like 
static dipoles. According to the arguments of sect. 5 and the formulae (|5.19|) the potential energy 
for the interaction is therefore proportional to 

ECa-pJOS-gJCa-VJOa-V)^ = 2«(a.0/3)(a.V)(/9.V)i (9.1) 



Recalling from ( p.24j ) that the spins for the individual Skyrmions are given by 

J (1) = Aq and J (2) = A/3 (9.2) 

and the isospins by 

lW=AJ (1) and I^=BJ {2 \ (9.3) 

we can write the potential energy for the long-range interaction between spinning triplets of dipoles 

as 

C x 2/<I (1) ■ I (2) ) ( J (1) • V) ( J (2) ■ V) 4- (9.4) 

R 
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where C is an unknown constant of proportionality. Remarkably this potential has exactly the 
same isospin and spin dependence as the zero-mass version of the one-pion exchange potential 
(OPEP) familiar from nuclear physics fIB |. In the usual treatment of the nuclear two-body problem 
in the Skyrme model one projects the static two-Skyrmion potential (|5.15|) onto tensor products 
of quantum states of free Skyrmions (which model nucleons) |12| . One then also reproduces the 
OPEP, which is regarded as one of the successes of the Skyrme model. Here we have offered 
an alternative way of extracting the OPEP from the Skyrme model which takes into account the 
Skyrmions' spin and isospin already at the classical level. One still needs to determine the constant 
C, however. This requires more thought. In particular one needs to calculate the dipole moments 
of a rapidly spinning Skyrmion. Centrifugal effects could change them significantly relative to 
their static values. 



10 Conclusion 

In this paper we have shown that a moduli space inspired technique (the relativised product ansatz) 
and a point-particle model (the dipole picture) lead to the same 12-dimensional Lagrangian system 
for asymptotic two-Skyrmion dynamics. The interpretation of the agreement between the two 
approaches is that well-separated Skyrmions interact like two point-particles with suitable scalar 
dipole moments (adjusted to fit the asymptotic field of a single Skyrmion) not only statically but 
also dynamically. 

Relativistic effects such as Lorentz contraction were important in our calculation, and to in- 
corporate them we had to modify the usual moduli space prescription. In fact, on the basis of our 
calculations one might be surprised that in soliton theories of Bogomol'nyi type the moduli space 
approximation without relativistic corrections satisfactorily reproduces the results of numerical 
calculations in the full field theory and agrees asymptotically with point-particle approximations. 
The reason for this, however, is a special property of soliton theories of Bogomol'nyi type: in such 
theories all fields in the moduli space are at the minimum of the potential, and therefore Lorentz 
contraction leads to corrections of order (velocity) 4 . By contrast, in theories where there are static 
inter-soliton forces, Lorentz contraction leads to corrections of order (velocity) 2 , i.e. of the same 
order as the kinetic energy. In such theories it is therefore not consistent to include the kinetic 
energy in a truncated Lagrangian but to ignore relativistic effects. 

While we could modify the moduli space approximation for two-Skyrmion dynamics, as pro- 
posed in 0, to take into account Lorentz-contraction of the individual Skyrmions, we have argued 
that it has to be abandoned altogether when studying the long-range forces between spinning 
Skyrmions. In its stead the dipole picture should be used to calculate those forces. Moduli space 
techniques are only reliable when the Skyrmions' separation and their angular speeds satisfy the 
inequality (|7.32|) or when the Skyrmions are so close together that their interaction forces them to 
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move collectively rather than individually. However, for near- coincident or coincident Skyrmions 
the product ansatz is no longer adequate and the more sophisticated methods of ref. || will be 
required. The Lagrangian L 2 would be useful in this context because it could interpolate between 
the moduli space and the point-particle approximation. 

The immediate question that arises from our work concerns the precise nature of spinning 
Skyrmions. Here we have argued that the asymptotic field is that of a spinning triplet of dipoles, 
but it would be interesting to study spinning Skyrmions numerically to check our argument and 
to calculate the field near the Skyrmion's centre. Such calculations should be carried out with 
massive pion fields, which is necessary for the existence of a stable rotating solution of the Skyrme 
equations. The results would flesh out or refute the argument of sect. 9, that it is the forces 
between spinning Skyrmions which one needs to understand in order to calculate nuclear forces 
from the Skyrme model. 
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